Study of the collapse of granular columns using DEM numerical simulation 
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Numerical simulations of the collapse and spreading of granular columns onto an horizontal plane 
using the Contact Dynamics method are presented. The final shape of the deposit seems to depend 
only on the aspect ratio a of the columns; these results are in good agreement with previous experi- 
mental work. In particular, the renormalised runout distance shows a power law dependence on the 
aspect ratio a, which is incompatible with a simple friction model. The dynamics of the collapse is 
shown to be mostly controlled by the free fall of the column. The energy dissipation at the base of 
the column can be described simply by a coefficient of restitution. Hence the energy available for 
the sideways flow is proportional to the initial potential energy Eq. The dissipation process within 
the flow is well approximated by basal friction, contrary to the behaviour of the runout distance. 
The mass ejected sideways is showned to play a determining role in the spreading process. As a 
increases, the same fraction of initial potential energy Eo drives more mass against friction. This 
additional dissipation give a possible explanation for power-law dependence of the runout distance 
on a. Beyond the frictional properties of the material, we show that the flow characteristics strongly 
depend on the early dynamics of the collapse. We propose a new scaling for the runout distance that 
matches the data well, is compatible with a friction model, and provide a qualitative explanation to 
the column collapse phenomenology. 



I. INTRODUCTION 

Although they have been a central subject in recent 
research, granular flows remain intriguing in many 
aspects of their behaviour [see0,0|. The understanding 
of such flows has obvious application in industrial 
processes, which often handle all kinds of powders and 
granules. But it is also crucial in the geophysical issue 
of catastrophic flows. Rock avalanches, landslides or 
pyroclastic flows all involve a granular solid phase, the 
influence of which can become momentous on the mean 
behaviour of the flow. 

In the absence of a clear physical background for the 
modelling of these large scale natural granular flows, 
their characterisation mainly relies on the observation 
of the flnal deposits, and in particular, of the final 
runout distance [see 0, 111- A simple effective basal 
friction fie is often advocated as the most convenient 
description of the dissipation process. This allows for 
a in-situ quantification of the mobility of the flow, 
knowing the runout distance and the initial height of the 
material. Moreover, basal friction is easily incorporated 
in continuous modelling, like shallow water equations for 
instance [see[3,l30|- However, the dependence of this 
effective friction on the nature of the material, on the 
characteristics of the collapse or on the dynamics of the 
flow, has yet to be specified. 

Systematic studies of the frictional properties and the 
spreading of granular flows on incline planes have been 
carried out [seeQ,^l|. However, until late, the simple 
case of a granular mass collapsing onto an horizontal 
plane had not been addressed. Recently, experimental 
works have tackled this problem, and studied the 



collapse and spreading of a suddenly released column of 
grains onto an horizontal plane [see[l|,|3,[i3. Essentially, 
the effect of the initial geometry of the column on the 
geometry of the final deposit has been investigated. 
The main result consists in scaling laws for the runout 
distance. In particular, when the initial aspect ratio 
of the collapsing column is high enough, the runout 
distance renormalised by the initial radius of the column 
shows a power law dependence on the initial aspect 
ratio. Moreover, this dependence varies depending on 
the conditions of the experiment, and in particular 
between the axisymmetric or quasi 2D configurations. 
Although simple models relying on Coulomb failure 
analysis or on shallow water approximation for the flow 
have been proposed, no clear physical understanding 
of the granular collapse process has been achieved yet. 
In particular, the scaling laws obtained for the runout 
distance are incompatible with a simple basal friction 
model, and question the mode of energy dissipation 
occurring in the successive steps of the collapse dynamics. 

In this context, the use of the numerical simulation to 
reproduce numerically the collapse of granular columns 
is expected to give new interesting insights in the prob- 
lem. Indeed, Discrete Element Methods allow for the 
simulation of each grain forming the granular mass, giv- 
ing access to each grain's trajectory. We can thus hope 
to access a picture of the collapse phenomena otherwise 
very difficult to obtain from experimental setup. 
The aim of the present work is to determine which are the 
mechanisms controlling the spreading dynamics: fall of 
the column, dissipation at the base, ejection of mass side- 
ways and dissipation within the fiow. Therefore, we have 
applied the Contact Dynamics algorithm [see[a,[l3| in two 
dimensions. The numerical procedures are explained in 
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sectional After briefly recalling and commentingthe ex- 
perimental results obtained by other authors [seeflLlMllOll 
in section Hnl the different regimes of spreading are quali- 
tatively described and scaling laws are established in sec- 
tion We observe a very good agreement between the 
experiments and the simulations. Details of the dynam- 
ics of the vertical fall and of the sideways spreading are 
shown in section The dissipation of energy, and the 
transfer of energy from vertical fall to horizontal motion, 
are examined in section IVTI We show that the amount 
of energy available for the spreading is simply propor- 
tional to the initial potential energy of the column and 
that basal friction is a very good approximation of the 
dissipation within the sideways flow. Eventually, we con- 
clude that the ability of the column to eject the grains 
sideways is a major aspect of the dynamics. This leads 
us to question the scaling laws obtained for the runout 
distance as perhaps fortuitous and corresponding to a 
transient regime. A new empirical fit, compatible with 
a friction law, and qualitatively describing the collapse 
phenomenology, is proposed in section fV ill Summary of 
the results and further discussion are presented in sec- 
tion EIIIl 



II. NUMERICAL PROCEDURES 

A. Simulation Method 

The numerical methods used for the simulation of 
granular material are known under the generic name 
of Discrete Element Methods. They take into account 
the individual existence of each discrete grain forming 
the media, and most usually neglect the role of the 
interstitial fluid filling the space between the grains. 
Such an approximation makes DEM methods mostly 
appropriate for the modelling of dry granular matter. In 
the absence of any dynamics induced by the surrounding 
fluid (the air in the present case), the behaviour of 
the collection of grains is entirely driven by the usual 
equations of motion, and the contact laws ruling the 
collisions between the grains. 

The modelling of the contact phenomena consists in 
a set of relations between the contact force and various 
quantities referring to the physical and/or chemical pro- 
cesses taking place at the surface of the two bodies in 
contact. These microscopic interactions can be of dif- 
ferent natures. They can involve for instance elastic or 
plastic deformation of asperities, adhesion due to Van der 
Waals forces, aging process... In any case, the complex- 
ity of these phenomena cannot be directly incorporated 
in the contact models, first due to practical reasons of 
numerical feasibility or efficacy, and then, because such 
complex models are nevertheless not assured to be re- 
alistic. Basically, DEM methods assume small elastic 
deformations, with possible viscous effects, and frictional 
dissipation as contact phenomenology 2]. 



It is beyond the scope of this paper to give a detailed 
account of the numerical method used, and further infor- 
mations will be found in dedicated work. Hence we will 
just specify the main hypothesis ruling the behaviour of 
the numerical grains. In the absence of a clear physical 
background for the modelling of contact phenomena, a 
possible strategy is to assume that grains are interacting 
through hard-core repulsion and non-smooth Coulombic 
friction only. This modelling choice is adopted in Con- 
tact Dynamics (CD) algorithm [see 0, 0| that we have 
applied for the present work. This implies that two grains 
have to touch for the contact force to be non-zero, and 
no distance interaction is permitted. Once a contact is 
formed between two grains, the latter cannot get closer, 
and any normal relative motion is repulsive: the grains 
are perfectly rigid. The microscopic coefficient of friction 
/i characterizes the Coulombic friction threshold ±nN, 
where N is the normal force at contact, and controls the 
friction dissipation process at grain scale. The tangential 
force T between two grains in contact can either be below 
the Coulombic friction threshold, and in that case no tan- 
gential slip motion is possible, or it can be exactly equal 
to the Coulombic friction threshold, and in that case, slip 
motion will dissipate energy. These contact laws allow for 
the modelling of multiple collisions between the grains, 
namely the fact that the grains are in contact with one 
or more neighbours while they are undergoing collisions 
with other grains. However, the simple case of a binary 
collision is also perfectly described. A Newtonian coef- 
ficient of restitution p is introduced, which controls the 
velocity of the grains after the collision knowing the ve- 
locity before the collision. This coefficient of restitution 
also appears in the prediction of the velocity of all the 
grains in the solving process. 

These contact laws are simplistic with regards to the mi- 
croscopic reality of contact phenomena. Nevertheless, 
they are sufficient ingredients to reproduce the collective 
dynamics of a collection of grains. 



B. Numerical Experiments 

Using the CD method, we have simulated two- 
dimensional collections of non-adhesive rigid circular 
grains. The diameter d of the grains is uniformly 
distributed in a small interval such as dmin/dmax = 2/3. 
This slight polydispersity of grains size is chieffy intro- 
duced to break any crystal like ordering of the grains 
which may have a non- negligible effect in 2D simulations. 
In the following, d denotes the mean grains diameter. 
The coefficient of restitution p at collision is the same 
between all the grains and between the grains and the 
bottom plane, and was chosen p = 0.5. As well, the 
microscopic coefficient of friction is the same at all 
contacts and equals p — 1. 

The numerical experiment consists in releasing a col- 
umn of grains in the gravity field onto a flat bottom plane. 



3 



FIG. 1: Preparation of a column of grains by a random rain 
of grains in the gravity field (top left picture). The column 
is characterized by its initial radius Ro and it initial height 
Hq. After the collapse, the final deposit is characterized by 
the runout distance Roo and the height Hoo 



and study the collapse and the resuhing spreading dy- 
namics of the granular mass. The bottom plane in our 
experiment is perfectly smooth. Laboratory experiments 
carried out on rough or smooth surfaces have shown that 
the general behaviour of the mass of grains (in particular 
the final runout distance) is not affected by the roughness 
of the plane, save a slight transformation of the deposit 
final shape [seeQ. 

The initial columns are prepared by means of a random 
rain of grains between two vertical walls. The compacity 
of the packing is cq ~ 0.82. The dimensions of the col- 
umn are its radius Rq and it height Hq, and a = Hq/Rq 
is the initial aspect ratio. At time t — 0, the vertical 
walls are instantaneously removed, and the column col- 
lapses due to gravity. The dimension of the final deposit 
are the final runout distance Roo, and the final height 
Hoo- The compacity of the final deposit is Coo — 0.78, 
namely close to the initial compacity in spite of a slight 
loosening of the packing of grains. This successive steps 
are illustrated in Figure 

We have carried out 25 simulations, with a ranging 
between 0.21 and 17, and counting from 1000 to 8000 
grains. 



distance. Using the notation mentioned above (namely 
Hq and Rq are the initial height and radius of the column 
respectively, a the initial aspect ratio, and Hoo and Roo 
are respectively the height of the final deposit and the 
runout distance), the first authors [see|l^ find: 
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Ro 



Rf) 



1.24 a, a<1.7 
1.6ai/2 a > 1.7. 



while Lajeunesse et al [seeQ obtain: 
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1. a, 



< 0.74 



2. a > 0.74. 



Moreover, quasi-2D experiments were carried out by 
Lube et al [see 1 1 If by releasing granular columns confined 
between two vertical walls. The following scalings were 
thus obtained: 



Ro 



Ro 



Ro 



1.2 a, a<2.3 
1.9a2/3 a > 2.3. 



Quasi-2D experiments were also carried out by Balm- 
forth & Kerswell [see ij , where the influence of the gap 
between the two vertical walls confining the column col- 
lapse is also addressed. The results indicates that the 
exponent of the power law is dependent on the size of 
the gap: 
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^^,0.65 Naj-j-Q-^^ gap 
Aa°-^ Large gap 



In these experiments moreover, the pref actor A varies 
depending on the material used, while previous authors 
found a mono-valuated prefactor, presumably due to 
the narrow range of experimented material in this latter 
cases. 

However, scalings found for quasi-2D experiments in a 
narrow gap co nfig uration give similar results for the two 
authors [seeQ,^, giving roughly Roo — RoRo oc 



,2/3 



The origin of the exponents is still open to discussion. 
No model yet has achieved a comprehensive explanation 
of the collapse dynamics. In particular, a simple friction 
model cannot account for them. Indeed, supposing that 
the initial potential energy of the column is completely 
dissipated by the work of the friction forces along the 
runout distance leads to: 



III. COMMENTS ON THE EXPERIMENTAL 
RESULTS 

Recently, experiments by Lube et al. [see'lO^ and La- 
jeunesse et al [see 8] have investigated the axisymmetric 
collapse of a column of grains onto a flat horizontal plane. 
The main result consists in scaling laws for the runout 
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where mg is the total mass of grains, and fj,^ is the effec- 
tive coefficient of friction (constant by definition). The 
existence of the exponents thus leads to the conclusion 




FIG. 2: Grains movement in squat columns. The snap- 
shots show successive instants t/Tac = 0, t/Taa = 0.25, 
t/Tac ~ 0.37 and t/Too = 1, where Too is the total dura- 
tion of the collapse. The column's aspect ratio is a = 0.37. 
In black are represented the grains the cumulated horizontal 
displacement of which exceeds the mean grains diameter d. 
The scale in the four pictures is the same. 




FIG. 3: Grains movement in tall columns. The snap- 
shots show successive instants t/Tao = 0.1, t/Taa = 0.25, 
t/Toc ~ 0.37 and t/Too ~ 1, whore Too is the total duration of 
the collapse. The column's aspect ratio is a = 9.1. In black 
are represented the grains the cumulated horizontal displace- 
ment of which exceeds the mean grains diameter d. The scale 
in the four pictures is the same, and the deposit has been 
truncated in the last picture. 



FIG. 4: Grains movement in intermediate size 
columns. The snapshots show successive instants i/Too = 0, 
t/Too = 0.25, t/Too = 0.37 and t/Too = 1, where Too is the 
total duration of the collapse. The column's aspect ratio is 
a = 0.9. In black are represented the grains the cumulated 
horizontal displacement of which exceeds the mean grains di- 
ameter d. The scale in the four pictures is the same. 




FIG. 5: The inner static cone. Final deposits of columns 
with a = 0.9, a = 0.73, a = 0.55 and a = 0.37. The inner gray 
cone coincides with the grains whose cumulated horizontal 
displacement is smaller than 5d, and shows a slope ~ 35° in 
all cases. The scale in the four pictures is the same. 



IV. SCALING LAWS FOR THE FINAL 
DEPOSIT 



that the dissipation process in the collapsing columns 
cannot be simply interpreted as a simple basal friction, 
and that the overall dynamics of the spreading might be 
more complex than usually postulated. 



A. Qualitative Description 

From a simple qualitative observation of the dynamics 
of the collapse, two different regimes can be distinguished 
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depending on the value of the initial aspect ratio a. 
In the first regime, for small a, the flow simply consists 
in the fall of the edges of the initial column. The motion 
is propagating from the edges inward, while a slope 
progressively builds up, along which the grains eventu- 
ally stabilize. In this regime, only the grains situated 
at the sides of the columns fall and flow as a result; by 
contrast, the grains situated inside the column have no 
motion and play no role at all in the spreading. This 
situation is illustrated in Fig |21 where four snapshots 
of a collapsing pile with a — 0.37 are displayed. In the 
course of time, we have represented in black the grains 
the cumulated horizontal displacement of which exceeds 
the mean grains diameter d. We observe that a majority 
of grains experiences smaller or no displacement at all, 
and that the upper surface remains undisturbed in its 
larger part. 

The second regime, namely for high a, is radically 
different. In that case, the whole column falls in a 
vertical motion in response to gravity, causing most of 
the grains to take part to the overall dynamics. Four 
snapshots of the collapse of a column are shown in FigO 
for a = 9.1. Again, are represented in black the grains 
whose cumulated horizontal displacement exceeds d. 
Only a small fraction of grains situated in the center of 
the column remain undisturbed. When a tends to oo, 
this fraction tends to zero. 
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FIG. 6; Final renormalised runout distance (-Roo — Rq)/Ro 
as a function of the columns aspect ratio a. 
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An intermediate case between these two regimes is 
shown in Fig^ for a = 0.9; the entire upper surface is 
affected by the sideways flow, but a well defined inner 
cone remains static. The shape of this inner cone is 
likely to be related to the frictional properties of the 
material, as suggested by 0. In Fig |S1 final deposits 
are represented for columns with a = 0.9, a = 0.73, 
a — 0.55 and a — 0.37. To make the distinction more 
obvious between the flowing region and the static or 
creeping inner part of the column, a stronger criterion is 
used to distinguish the flowing grains: their cumulated 
horizontal displacement must exceed 5d. By this means, 
the slope of the inner cone becomes very neat, and allows 
for the comparison of the deposits. We observe that the 
slope of the inner cone remains nearly constant, namely 
around 35°, independently of a. It is very tempting 
to assume that this slope reflects the packing mean 
frictional properties rather than the flow dynamics, and 
that the collapse, at least for small values of a, results 
from a Coulomb failure as assumed by Lajeunesse in 
In this case, assuming an hydrostatic stress gives us 
an internal angle of friction ip c± 20° for the packing 
of grains, equivalent to a coefficient of friction 0.36. 
We note that this value does not compare with the 
coefficient of friction acting at contacts between grains 
/J, = 1, and neither with the effective coefficient of basal 
friction later evaluated in section IVTI For greater values 
of a, the static inner cone is destroyed by the vertical 
dynamics of the upper grains and can no longer be 



FIG. 7: Final maximum height of the deposit renormalised 
by the column initial radius Hac/Ro (left), and mean height 
of the final deposit renormalised by the column initial height 
Hoo/Ho (right) as a function of the columns aspect ratio a. 



observed in the final deposit. 

The general observations on the shape of the col- 
lapsing columns just exposed are very close to previous 
experimental description of the collapse phenomenology 



B. Scaling Laws for the Final Deposit 

The shape of the final deposit resulting from the col- 
lapse and the spreading of the granular column is first 
characterized by the final runout distance i?oo • For each 
experiment, Roo is evaluated from the position of the 
grains connected to the main mass by at least one con- 
tact. In other words, grains ejected from the flow, and 
undergoing a solitary trajectory independently of the col- 
lective behaviour of the flow, will not be taken into ac- 
count. The final height of the central conical region Hoc 
is also evaluated; it corresponds to the highest point of 
the deposit. Moreover, we evaluate the mean height of 
the deposit Hoc from the total area covered by the grains. 
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This last quantity is equivalent to a measure of the final 
potential energy of the deposit. 

The evolution of the distance run by the grains nor- 
malised by the initial radius of the column (i?oo ~Ro)/ Ra 
is plotted in Fig|Slas a function of a. We observe the fol- 
lowing dependence: 



Ro 



Rq 



Ro 



2.5 a, 
3.25 a° 



<2 
> 2. 



As observed in laboratory experiments, the runout dis- 
tance obeys two different behaviours depending on the 
value of a. For small values of a a linear dependence is ob- 
served, while for larger a, the dependence is a power law. 
The scalings we obtain are in very good agreement with 
the scalings observed from quasi 2D laboratory experi- 
ments, for which the exponent observed for large enough 
a IS 2/3 [seed lllj . and the transition between the two 
behaviours occurs at a = 2.3 [see lllj . The prefactors 
2.5 and 3.25 obtained in the numerical experiments are 
higher than those observed experimentally by Lube et al, 
namely 1.2 and 1.9. This difference is very likely due 
to the respective frictional properties of the material; the 
circular shape of the numerical grains is for instance a fac- 
tor of mobility. Moreover, the dissipation induced by the 
friction with the two vertical walls confining the grains 
in quasi-2D configuration is certain to affect the value of 
the prefactors. 

The normalised final height Hoa/Ro is plotted against a 
in Fig El we observe: 



Hoc 

Rq 



1.6aO-35, 
1.45, 



a < 10 
a > 10. 



The 2D laboratory experiments show a similar behaviour 
(with an exponent 0.4 [see^| or 0.5 [seeQ), but no tran- 
sition is observed for a ~ 10. 

The increase of the quantity of grains spreading sideways 
appears clearly when plotting the mean height of the de- 
posit normalized by the initial height Hoof Hq against the 
initial aspect ratio a. From Fig [3 the dependance shows 



1, a<0.8 
0.4 a^" ^ a>0.i 



as can be expected from the scaling of the runout dis- 
tance and the mass conservation RqoHoc ~ RqHq. This 
relation shows the increasing transfer of potential energy 
to sideways spreading motion, and the decrease with a of 
the potential energy of the final deposit. 
As discussed in the introduction, these scaling laws are 
incompatible with a simple frictional behaviour which 
would bring Rqc/Rq oc a and Hoo/Rq oc a^^. A first 
hypothesis is that the dynamics of the grains at the bot- 
tom of the column is responsible for a complex dissipa- 
tion process dependent on the initial aspect ratio a. As a 
consequence the energy available for the spreading would 
also depend on a, and possibly cause the dependance of 
{Roo — Ro)/Ro with a to be a power law. This aspect 
will be discussed further in section IVTl 
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FIG. 8: Dependance of (tf/to)^ on the column initial aspect 
ratio a, where t/ is the time of free fall of the top of the 
column and to — {2Ro/g)^^'^ ■ 



V. DYNAMICS OF THE COLLAPSE AND 
SPREADING 

A. The Vertical Fall 

The dynamics of the collapse is first induced by the 
vertical fall of the grains. However, as stressed in sec- 
tion llV K\ different behaviour can be observed depending 
on the value of a. Computing the position h of the top of 
the column in the course of time, we are able to compare 
h with the free fall position given by Hq — 0.5gt^, and to 
evaluate the time t f during which the top of the column 
is in free fall. The criterion used for the free fall to cease 
is \h — {Hq — 0.5gt'^)\ < d (where d is the mean grains 
diameter) . 

In FiglHlwe have plotted {tf/to}'^ as a function of the as- 
pect ratio a, where to = {2Rq/ gY/^ . For small values of 
a, the time of free fall tj is nearly zero. However, when 
a increases, the following relation is satisfied: 



to 



0.95a- 2.5, 



/2(i7o-2.5i?o) 
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This means that the top of the column is undergoing 
a free fall along a distance of nearly iJp — 2.5Rq. In 
other words, columns with a ^ 2.5 are subjected to free 
fall, while columns with a ^ 2.5 are not. We believe 
this transition in the vertical dynamics of the column 
at a ~ 2.5 to be at the origin of the transition observed 
in the scahng law relating (i?oo ~ Ro)/Ro to a (see 
FigEJ. 
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FIG. 9: Simultaneous collapse for a — 8.8 and a — 17.4. 
The eight top pictures show the indifference of the spreading 
on the height of the columns during the free fall dynamics. 
The eight bottom pictures show how the additional material 
(in black) pushes aside and eventually covers the underlying 
material (in gray). 



The fact that the top of the column is in free fall 
entails that the upper part of the column is not affected 
by the complex spreading process going on at the 
bottom. As a consequence, two columns with the same 
initial radius Rq, but two different initial heights Hq and 
Hq, should behave the same as long as the top of the 



smallest column remains above « 2.5i?o- The spreading 
at the base should not be affected by the height of the 
column above, and the top of the column should no 
see the spreading process underneath. Once this limit 
height « 2.5Rq is reached, the smallest column will soon 
finish its course, while the higher column will be further 
accelerated. 

This behaviour is clearly visible in the series of pictures 
shown is Fig representing the simultaneous collapse 
of two columns with heights Hq and Hq ~ 2Hq respec- 
tively, and same initial radius Rq. In the second column, 
the grains initially situated above the height Hq are rep- 
resented in black to allow for the comparison of the two 
dynamics. From the beginning until Hq ~ 2.5i?o (first 
four pictures from top to bottom), the top of the two 
columns remain undisturbed, and the spreading process 
occurring at their base is exactly identical. Then, we can 
observe, in the respective evolution of the mass of grains 
initially situated under the height Hq (represented in 
gray), the effect of the fall of additional grains (in black) 
over the underlying deposit. In particular, the dynamics 
of pushing aside grains that would otherwise remain 
in the vicinity of the bottom of the column is obvious. 
Eventually, the black grains cover the underlying gray 
ones. 



B. The Sideways propagation 

In the course of time, the sideways flow propagates 
outward; we denote r the front position at any time t; 
eventually, r reaches the final value Roo- The total du- 
ration of the collapse at which the sideways propagation 
stops and the whole deposit comes to repose is denoted 
T 

In order to compare the dynamics of the spreading 
for different values of a, we plot in the same graph in 
Figure ^| the evolution of the front position normalised 
by the final runout distance (r — Ro)/Roo as a function 
of the time normalised by the total duration of the 
propagation t/Too, for a = 0.9, a = 1.8, a = 3.1, a = 7.7, 
a = 10.7 and a = 15.7. The plots nicely collapse into a 
master curve, showing first a period of acceleration of 
the front, followed by a constant propagation regime, 
and then a deceleration period. 

From the initial geometry of the column two character- 
istic times to = (2i?o/g)^/^ and Tq = {2Ho/gY/'^ can be 
formed, corresponding to the time of free fall along the 
distance i?o and H^ respectively. The collapse duration 
Too is plotted against Tq in graph^] We observe a linear 
relation 

f2H \ 

Too ^ 2.25To = 2.25 f — ^ j , 

implying that the duration of the experiment is con- 
trolled by the free fall of the column of initial height 
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FIG. 10: Front position renormalized by the final runout 
distance (r — Ro)/Roo as a function of the time renormahsed 
by the total duration of the collapse t/Too for different values 
of a. 




FIG. 11: Total duration of the collapse Too as a function of 
To = {2gHoY''^. The linear relation gives Too ^ 2.25To. 



Hq. This is in agreement with experimental results, for 
which the duration of the experiments is found to be 
nearly STq [see^. 

The period of acceleration of the propagation front is 
well characterized when plotting for different values of a 
the normalised front position (r — Rq) / Rq as a function 
of the renormaliscd time t/to ('Figll2|). For high values of 
a, we clearly distinguish the acceleration phase followed 
by a constant velocity propagation phase. This evolution 
is not so obvious for small values of a, for which the de- 
celeration phase occurs early and leaves less time for a 
constant velocity regime to settle. 

When plotting (r — Ro)/Ro against a in a log-log repre- 
sentation, we see, up to t/tg — 1.5, and for a 1.8, that 
the following relation is acceptable: 



^ ^ 0.68 ( - 

Ro \to 



Although this approximation is made over a very short 
time interval, it suggests that the onset of the spreading 
is driven by the free fall dynamics. 

Where a constant velocity regime can be observed, the 
following relation is satisfied: 



r-Ro 



to 



which is equivalent to 



r ~ vot, r > 2i?o- 

where vq — (2gi?o)^^^ is the front propagation velocity, 
once a constant velocity regime is settled, after the front 
position has already run a distance 2Rq. 

Since the column is undergoing a free fall for a ^ 2, 
the cumulated mass of grains m(t) expelled as a result 
of the collapse is given by m(t) cx psRogt^, where ps is 
the surfacic mass of the grains packing. As the greater 
part of the front propagation involves a constant veloc- 
ity, r ~ VQt, the increasing mass debit at the bottom of 
the column can only be accommodated by an increase 
of the height of the sideways flow. Moreover, the grains 
reaching the bottom, after they have been accelerated 
in the gravity field, have a greater momentum than the 
grains preceeding them. This effect is responsible for the 
existence of a wave propagating outwards and transfer- 
ring the mass from the center towards the margin of the 
spreading for high values of a. An extreme case of this 
"mass propagation" phenomena is illustrated in Fig 1141 
where the sideways flow is represented for a = 70. 

This effect is more important in 2D configurations than 
in axisymmetric ones, for which the increase of the sur- 
face available for the spreading is quadratic with the 
front propagation, while the front propagation has been 
shown to obey the same behaviour r — vqI [see,^. This 
suggests a purely geometrical explanation of the differ- 
ence observed in the scaling laws in axisymmetric and 
2D experiments. Indeed, the front propagation veloc- 
ity vq — y/{2gRQ), and the typical time of the experi- 
ment To — ^/{2Ho/g), give for the runout distance the 
straightforward scaling law: 



{Roo — Ro) — voTo — {Hf)Rf) 



,1/2 



corresponding to the observation (Rao ~ Ro)/Ro — a^^^ 
in axisymmetric collapses. The difference with the 2D 
case (namely {Rod — Ro)/Ro — a^^^) could be due to 
the larger increase of mass with important momentum 
in the sideways flow, whose effect would be to lengthen 
the deceleration phase. The contribution of the decel- 
eration phase to the runout distance is non-negligible in 
2D. Three examples of the contribution of the decelera- 
tion phase to the final runout distance are displayed in 
Figure El up to one third of the total runout distance is 
gained during the deceleration. This effect is presumably 
more important than in axisymmetric configuration, for 
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FIG. 13: Front position renormalised by the column initial 
radius (r — Ro)/Ro as a function of the renormalised time 
t/to for different values of a. The linear evolution is stressed 
by the dashed line. Dashed-dotted lines indicate roughly the 
period of deceleration for each curve. 



which the mass debit is more easily accommodated by 
the expanding surface of the flow. This could explain the 
difference of the exponents (1/2 in axisymmetric and 2/3 
in 2D) in the scahng laws. A systematic analysis of the 
deceleration phase should help to solve this issue; it is 
however not undertaken in the present paper. 



FIG. 14: Successive snapshots of the spreading of a column 
with a very high aspect ratio a = 70. We observe the transfer 
of mass from the center of the collapse towards the front of 
the flow. The scale is the same on all pictures. 



VI. ENERGY TRANSFER AND DISSIPATION 

A. Time evolution 

Basically, three successive stages can be identified in 
the history of a grain falling within the column. In a first 
step, its initial potential energy is converted into vertical 
motion, and if the initial height of the grain allows for 
it, it will be accelerated down to the bottom. There, in 
a second step, the grain will undergo collisions with the 
bottom plane or the surrounding grains, and its vertical 
motion will be converted into horizontal motion. In a 
third step, the grain eventually leaves the column base 
area and flow sideways. Of course, this process involves 
a collective dynamics of collisions and momentum lost 
and transfer, whose complexity makes any prediction 
of the trajectory of any grain dubious. For the same 
reason, a high initial potential energy is no individual 
guaranty, for a single grain, of a long sideways travel. 
As an illustration of the complexity of the vertical to 
horizontal motion transfer, successive snapshots of the 
deformation of a collapsing column with a ~ 9.1 are 
displayed in Fig^l The grains initially situated at the 
margins of the column, in the central area, and at the 
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top, are represented in black. In the course of time 
we observe that the grains traveling farther are not 
those initially at the top. On the contrary, grains which 
were at middle height or even lower are nearer to the 
spreading front area. This behaviour is also visible in 
the deformation of the inner black slice. 

The conversion of momentum from vertical to hori- 
zontal direction is likely to depend on the value of a. In 
Fig 1161 we have represented for two collapsing columns 
with a = 2.6 and a = 15.7, the time evolution of the 
potential energy Ep, the vertical kinetic energy and 
the horizontal kinetic energy E^^ normalised by the ini- 
tial potential energy Eq, with 



E„ 



E., 



p=i 
1 



where i — {x,?/}, Np is the total number of grains, nrip 
their mass, hp their height and Vp their velocity. 
From these graphs we first see that a higher proportion 
of the initial energy is dissipated by the column with 
a = 15.7, in agreement with the scaling of Hoo/Hq in 
Figure |3 which is equivalent to the scaling of the final 
potential energy to the initial potential energy. Then, we 
observe that the conversion of energy from potential to 
kinetic in the vertical direction is much more efficient for 
a = 15.7. By contrast, a greater proportion of vertical 
kinetic energy is transferred into the horizontal direction 
for a = 2.6. This suggests that the ability of the column 
to use its initial energy for spreading might depend on a. 



B. Conversion from Potential to Kinetic Energy 

For each value of a, we compute the kinetic energy 
of the system averaged over the total duration of the 
collapse Too in the vertical and horizontal direction: 



J- DO 



p=l 



where i = {x,y}. To compare these values with pre- 
dictions of the mean kinetic energy over the free fall 
time, and using the fact that Too — 2.5ro, with Tq = 
y/{2Ho/g), the mean kinetic energies {Ek^} and {Ek^) 
are multiplied by a factor 2.25 in the following discus- 
sion. 

The mean kinetic energy of a column of height Hq and 
radius Rq submitted to free fall over a distance Hq — fci?o, 
where A: is a constant, and averaged over the free fall time 
Tf = ^{2{Ha - kRa)/g) is given by: 

f-Tf 
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1 



/ Jo 





FIG. 15: Deformation of a collapsing column with a — 9.1, at 
t/Too = 0, t/Tac = 0.2, t/T^ = 0.28, t/T^ = 0.35, t/T^ = 
0.45, and t/T^ = 0.75. 
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FIG. 16: Time evolution of the potential energy Ep (plain 
line), the vertical kinetic energy Eky (dotted line) and the 
horizontal kinetic energy Ek^ (dashed line) of the column for 
a = 2.6 (top) and a = 15.7 (bottom). 
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if we assume no influence of the grains which have al- 
ready reached the bottom on the falling material. This 
eventually leads to the relation 



Eo ^ a'^l5 5 a' 



(3) 



For high values of the aspect ratio a (or for fc = 0), the 
mean vertical kinetic energy is bounded by 4/15. For 
intermediate values of a, {Eky) is expected to obey the 
relation|21with k = 2.5. The mean vertical kinetic energy 
{Ek ) renormalised by the initial energy Eq is plotted in 
Figfni as a function of a. In the same graph is plotted 
the prediction g[a) = (1 - 2.5/a)(4/15 + 2/5 x 2.5/a) 
from the relation O We observe that the simulations 
points and the prediction are very unlike, presumably 
because of the collective dynamics of the grains reaching 
the bottom of the column. However, following the pre- 
diction, the evolution of the mean vertical kinetic energy 
shows an increase (though much slower than would be 
expected) towards the limit value 4/15, which appears in 
the following approximation of {E^y ) : 



En 



15 



(l-a-"-"22). 



We define the coefficient of restitution py characteris- 
ing the transfer of potential energy into vertical motion, 
so that 

(Eky) = PyEa, 

where py — 4/15(1 — a^""^^) is an increasing function of 
a. For a cxd, py — > 4/15, namely the potential energy 
is converted into vertical motion following the ideal free 
fall case. For smaller values of a however, there exists a 
transient regime during which the vertical kinetic energy 
increases with a. Similar transient behaviour in the 
dynamics of the sideways flow is discussed in next section. 

The horizontal kinetic energy (Ek^), normalised by £^01 
is plotted in Fig ^1 as a function of a. We observe 



Eo 



0.16, 



for a > 2.5. As preceedently we define the coefficient 
of restitution p^ characterising the transfer of potential 
energy into horizontal motion: 

{Ek^) = PxEq, 

with px — 0.16. Surprisingly, in spite of the complex- 
ity of the process taking place at the bottom of the col- 
umn, the energy available for the horizontal motion is 
simply proportional to the initial potential energy, and 
only implies the introduction of a constant effective co- 
efficient of restitution. Hence we can no longer sus- 
pect the dissipation process at the bottom of the col- 
umn to be at the origin of the divergence with the sim- 
ple friction model described in equation [3 We have 
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FIG. 17: Vertical kinetic energy {Eky} averaged over the 
total duration of the collapse Too and normalised by Eo, as 
a function of the columns aspect ratio a. 
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FIG. 18: Horizontal kinetic energy (Ek^) averaged over the 
total duration of the collapse Too and normalised by aEo, as 
a function of the columns aspect ratio a. 



seen in section IIIII that the scaling laws implied that 
the relation mogHo = peinogRoo is not satisfied. Since 
{EkJ = PxEq, we cannot have {EkJ = HemQgRoa either, 
which would lead to Rao = 0.16//ieffo- We may thus sup- 
pose that the assumption of a basal friction controlling 
the energy dissipation in the sideways flow is wrong. This 
aspect is tackled in the next section. 
Computing Eky and E^^ the maximum kinetic 
energy in the vertical and horizontal direction during the 
collapse, the evolution of these two quantities as a func- 
tion of a is similar to the evolution of the corresponding 
mean quantities observed in Figure lTTI and lT^ The max- 
imum values are about three times the mean values, as 
could be expected from the time evolution of Eky and 
Ek^ shown in Figure El 



Characterizing the ability of the columns to convert 
vertical kinetic energy into horizontal kinetic energy by 
the ratio Px/ Py, we obtain a decreasing function of the 
aspect ratio tending towards px , suggesting again the ex- 
istence of a transient regime in the dynamics of spreading. 
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C. Basal Friction 

The basal friction is the dissipation mechanism most 
often postulated for the dense flow of model granular me- 
dia . an d is indeed proven to be a relevant description [see 
IliL Il5l | . It implies the definition of an effective coefficient 
of friction fi^, which is a mean phcnomcnological transla- 
tion of the more complex dissipation process taking place 
at smaller scale through collisions and contact friction. 
To check whether basal friction is a satisfactory descrip- 
tion of the dissipation occurring within the sideways flow, 
we compare the exact amount of energy made available 
for the flow with the work done by the mass of grains over 
the distance they eventually run. Therefore, we define 
in the collapse two vertical sections S situated at —Rq 
and Rq, and delimiting the initial position of the mass 
of grains and the actual horizontal flow, as showned in 
FigEl Integrating over the total duration of the collapse 
Too , we compute the exact amount of energy Es crossing 
the sections S, and thus taking part to the flow, and the 
exact mass of grains mg going through S. The energy 
Es is computed from both potential and kinetic energies 
of the grains. In the case of a simple friction process in 
the sideways flow, we should observe, independently of 
the dynamics occurring within the interval [—Ro,Ro], 

Es = fiemsg{Rs - Ro), 

where Rs is the final distance run by the mass of grains 
nis, and is taken equal to the final position of the center 
of mass of ms- 

The energy Es normalised by the work of the mass 
ms over a distance Rq is plotted as a function of the adi- 
mensional flowing distance {Rs — Ro)/Ro in Figure EUk 
for all the collapse experiments. We obtain a neat linear 
dependance establishing that energy dissipation in the 
flow is very well approximated by basal friction. The 
value of the effective coefficient of friction /ie is given 
by the slope of the linear relation and is found to be 
Me ^ 0.47. 

If we assume that all the grains travel the final runout 
distance i?oo — Ro, we again obtain a linear dependence, 
as seen in Figure I^Ub. but with a much smaller effective 
coefficient of friction /Xg ~ 0.16 expressing the maximum 
mobility of the flow. In any case, the following relation 
is verified: 

Es = ^Jie'msg{Roo - Ra)- (4) 

In Figure 121b . the plot of Es /Eq a.s a function of a 
shows that as soon as a > 3, Es/Eq ~ 0.44. As could be 
inferred from the evolution of the mean horizontal energy 
(previous section), the energy available for the flow is 
simply proportional to the initial potential energy. The 
relation 0] can be rewritten: 

O.UEo ~ Hemsg{Roo - Ro), (5) 



or equivalently, 

O.MmogHo ~ fiemsg{Roc - Ro)- (6) 

Finally, this leads to the following expression for a > 3: 

{Rao - Ro) ms 
Rq mo 

This relation suggests that the disagreement between the 
scalings observed experimentally and the simple friction 
model (see expression might rest in the definition of 
the mass of grains actually flowing. 

The renormalised mass of grains actually crossing the 
section S and taking part to the flow ms/mo is plotted 
in Figure I^Tb as a function of a. There exist a function 
/ of the aspect ratio such as 

~ 1 - / O , 

Too 

with f{a) when a — > oo. No argument for the empiri- 
cal fit f{a) = 4.8(6-|-a)~^ is proposed; the main point of it 
is that it captures a major aspect of the flow phenomenol- 
ogy: the ejection of grains sideways. For small values of 
a, a small fraction of grains flows sideways, and most of 
the mass remains trapped in the interval [—Rq,Ro\, at 
the bottom of the initial column. The fraction of flow- 
ing grains increases with a. When a becomes high, this 
fraction tends towards 1, as can be seen in Fig 1141 where 
the mass flowing sideways in represented for a = 70. The 
increase of the proportion of mass actually taking part to 
the flow can be related to the increase of vertical kinetic 
energy relatively to the initial energy (see Fig [T7|l . 
The increase of mass flowing sideways causes the friction 
dissipation process to be more efficient; while a increases, 
the same fraction of initial potential energy Eq is driving 
more mass against friction. This additional dissipation 
may well explain why an exponent lower than 1 appears 
in the scaling law. 

This progressive increase of the fiowing mass, bounded 
by the initial mass toq, also suggests that the dynamics 
of the collapse, in the range of a experimentally investi- 
gated, is transient. 

VII. NEW SCALINGS FOR A TRANSIENT 
REGIME 

In the range of aspect ratios a investigated in this 
present work, as well as in preceedent experimental 
work [see Q, |^ ^3 ; the collapses of the columns mainly 
differ in the mass of grains ejected sideways. The higher 
the value of a, the greater the proportion of grains 
ejected. Because of this difference, the propagating flows 
do not involve the same proportion of the initial mass of 
grains. For increasing values of a, the increase of mass 
will cause the work done by the flow to be more efficient, 
affecting thus the runout distance. 
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FIG. 19: During the collapse, the exact amount of energy Es 
and the exact mass of grains mg crossing the section situated 
at — -Ro and _Ro are evaluated. The center of mass of the 
sideways flow is in Rs- 
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FIG. 20: Energy available for the flow Es renormalised by 
the work of the flowing mass ms over _Ro as a function of 
the renormalised distance run by the center of mass of the 
flowing grains (-Rs — B«)/Ro (graph a), and as a function of 
the renormalised runout distance of the flow (i?oo — Ro)/Ro 
(graph b) 



However, the mass of grains flowing tends towards the 
initial mass mo when a increases, so we anticipate the 
dependance of the dissipation process on a to vanish. In 
this limit, the amount of energy dissipated by the flow 
would only be dependent on the runout distance. 
The dependance of the proportion of mass actually flow- 
ing ms/mo on a shown in Figure [?Tb reflects this be- 
haviour and the asymptotic evolution towards a new 
regime. However the scaling laws displayed in Figure |^ 
and Figure [7| do not at all express the existence of this 
transition. 

On the basis of our observations of the phenomenology of 
the sideways flow, we can thus suggest that there exists 
a function /(a), verifying /(a) when a — > oo, such 
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FIG. 21: Energy available for the flow Es renormalised by 
the initial potential energy Eq as a function of a (graph a), 
and mass of grains ms actually taking part to the flow renor- 
malised by the initial the total mass of grain mo as a function 
of a (graph b). 



as: 



Roo — Ro 
Ro 



(!-/(«)) 



(7) 



We have represented in Figure [53 the evolution of 
(^oo ~ Ro)/Ro with a, the power law approximation 
3.25a'^'^, and the empirical fit a/(l — /(a)), where 
/(a) = 4.8(6 -I- a)~^ describes the grains ejection process 
(see section above). We observe that in the range of 
aspect ratios investigated, the last choice is as acceptable 
as a power-law dependence. 

No argument on the form of the function /(a) will be 
discussed at this stage. Basically, it represents the ad- 
ditional dissipation entailed by the increase of the pro- 
portion of mass flowing, while the proportion of energy 
available for the flow is constant (at least for a > 3 in our 
simulations). The choice of an approximation of the form 
a/(l — /(a)) for the runout distance has two implications: 

1. For large a, the runout distance will eventually in- 
crease like the height of the column, as expected in 
a simple friction dynamics. 

2. It suggests that the power-law dependance might 
be fortuitous. 

In the absence of a comprehensive model explaining 
either one or the other approximation, the proposition 
discussed here remains purely speculative. However, it 
matches nicely the numerical results, and provides a qual- 
itative explanation for the non linear behaviour of the 
runout distance with a. Finally, it suggests that the key 
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FIG. 22: Final renormalised runout distance (iioo — Ro)/Ro 
as a function of a. Two different approximations are plotted: 
the power-law fit in dashed line, and the function a/(l — /(a)) 
in dotted line. 



point of the collapse problem lies in the dynamics of ejec- 
tion of the mass from the initial column itself, rather than 
in the characteristics of the sideways flow. 



VIII. SUMMARY AND CONCLUSION 

We have numerically investigated the collapse and the 
spreading of two-dimensional columns of grains onto a 
vertical plane using the Contact Dynamics method. This 
approach allows for a detailed analysis of the dynamics 
of the collapse taking into account the individual grains 
energy and trajectory. Our results are generally in good 
agreement with previous experimental work carried out 
in quasi-2D configurations [see HI ITll| . 
The collapse is first described in terms of the shape 
of the final deposit, and more specifically in terms 
of runout distance. A power law dependence of the 
renormalized runout distance with the initial aspect 
ratio of the columns is found for high aspect ratios, and 
a linear dependence is found for low aspect ratios. These 
scalings, experimentally observed by previous authors, 
are incompatible with a simple friction model of the 
collapse dynamics. 

We show that the collapse is driven by the free fall of the 
column for high enough aspect ratios. The existence, or 
absence, of free fall dynamics, can explain the existence 
of two different scaling laws for the runout distance 
depending on the aspect ratio. The propagation of 
the front involves a constant velocity phase, followed 
by a deceleration phase of significant contribution to 
the runout distance. The analysis of the mean kinetic 
energy of the grains shows that the dissipation process 



occurring at the bottom of the column can be simply 
described by a constant coefficient of restitution. In 
particular, the energy avalaible for the sideways flow 
is simply proportional to the initial potential energy. 
The detailed analysis of the energy dissipated in the 
sideways flow and the work of the flowing mass clearly 
establishes that friction is a good approximation of the 
dissipation process. Finally, we point at the dynamics 
of mass ejection sideways during the column collapse as 
playing a predominant role in the spreading dynamics, 
and as being responsible for the non-linear behaviour 
of the renormalised runout distance. This allows us to 
suggest that the scaling laws preceedently discussed for 
the runout distance are fortuitous, and should no longer 
apply when the aspect ratio increases. A new empirical 
flt is proposed, which is compatible with a friction model. 

This conclusion has different implications: 

1. Mass ejection sideways is a mechanism strongly de- 
pendent on the geometry of the collapse. In par- 
ticular, we expect its effects to be more important 
in a 2D configuration, than in the case of an ax- 
isymmetric collapse. This might be at the origin 
of the difference is the scaling laws for the runout 
distance observed in 2D (or quasi-2D) and axisym- 
metric configurations. In the limit of high aspect 
ratios however, the mass ejection tends towards the 
ejection of the totality of the initial mass. In that 
limit, differences should no longer be observed be- 
tween the 2D and the axisymmetric configuration. 
In any case, our results suggest that an analytical 
expression of the runout distance should account 
for the process of the ejection of grains at the bot- 
tom of the collapsing column. 

2. The runout distance appears to be strongly depen- 
dent on the fall dynamics and not only on the effec- 
tive flow properties, namely effective basal friction 
in our case. Although high aspect ratios are diffi- 
cult to find in natural context, many rock falls or 
slope destabilisations involve a strong acceleration 
(and possibly free fall) , which is a key aspect of the 
material ejection. In a geophysical perspective, this 
suggests that the mobility of a natural flow, usu- 
ally defined as the ratio of the runout distance to 
the initial height of the material, is related to the 
early dynamics of the mass release as well as to the 
flowing properties of the material. 

3. Since the sideways flow undergoes a simple basal 
friction dissipation process, its modelling using 
shallow water approaches is possible, but not 
straightforward. Indeed, a difficulty lies in the de- 
scription of the initial conditions represented by the 
vertical column collapse, which escapes, intrinsi- 
cally, the shallow water assumptions. The issue 
is to achieve a correct description of the mass flux 
with a limitation of the energy induced by the col- 
umn fall. The column collapse dynamics, until now 
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correctly reproduced only for low aspect ratios [see 
0, , could also be recovered for high aspect ra- 
tios [seeQ. 

The influence of the material properties (inter-grains 
friction /i and restitution at collision p) on the overall 



dynamics of the collapse and the spreading will be the 
subject of further works. 

This work was supported by the Marie Curie European 
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